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We use Landau's theory of a normal Fermi liquid to derive expressions for the static response 
of a system with a general tensor interaction that conserves the total spin and the total angular 
momentum of the quasiparticle-quasihole pair. The magnetic susceptibility is calculated in detail, 
with the inclusion of the center of mass tensor and cross vector terms in addition to the exchange 
tensor one. We also introduce a new parametrization of the tensor Landau parameters which 
significantly reduces the importance of high angular harmonic contributions. For nuclear matter 
and neutron matter we find that the two most important effects of the tensor interaction are to give 
a contribution from multipair states and to renormalize the magnetic moments. Response to a weak 
probe may be calculated using similar methods, replacing the magnetic moments with the matrix 
elements of the weak charges. 



I. INTRODUCTION 

For systems with central interactions, Landau's the- 
ory of normal Fermi liquids provides an economical way 
of characterizing many low-temperature properties. The 
theory was applied to atomic nuclei by Migdal and collab- 
orators, and in that work it was generally assumed that 
the non-central contributions to the effective nucleon- 
nucleon interaction were small jj. The generalization 
of Landau theory to include effects of the tensor force 
was made by Dabrowski and Haensel OJ FJl Subse- 
quently, estimates of tensor contributions to the effective 
interaction were made in Refs. [5J, |y, |jj, |8| • 

The stimulus for the present work arose in the con- 
text of astrophysics. In the physics of collapse and the 
subsequent evolution of a neutron star, the properties of 
neutrinos in dense matter are a ke y in gredient M . As 
demonstrated in Refs. [13, EH E3, 113, H3, EE EH E3 . 
neutrino scattering and absorption rates are sensitive 
to nuclcon-nuclcon interactions, especially their spin- 
dependent parts. Direct calculation of effective interac- 
tions is difficult, but for systems with only central in- 
teractions, information about the effective interaction at 
long wavelengths may be obtained directly from a knowl- 
edge of static properties of the matter. For example, 
the spin susceptibility of a single-component Fermi liq- 
uid with two spin states is given by 

where fio is the magnetic moment of a particle in free 
space, N(0) = m*p F /n 2 h 3 is the density of states per 
unit volume at the Fermi surface, m* is the quasiparti- 
cle effective mass, pp is the Fermi momentum and Go 
is the Landau parameter describing the isotropic part 



of the spin dependent contribution to the quasiparticle 
interaction. Since calculations of the magnetic suscepti- 
bility of neutron matter exist Eg, ^ is relevant to ask 
to what extent it is possible to deduce properties of the 
effective interaction from such data. A related question 
is the extent to which tensor contributions to quasipar- 
ticle energies and interactions alter neutrino scattering 
rates, which were previously calculated neglecting tensor 
effects in Ref. 11]. 

As shown in earlier work E3> the tensor force influ- 
ences the static response of a system in a number of 
different ways. One is that the magnetic moment of a 
quasiparticle is different from its value for a particle in 
vacuum. In particular, the magnetic moment is not a 
scalar, as it is for systems with central forces only. As 
long ago as 1951, Miyazawa [2(| calculated explicitly how 
the magnetic moment of a nucleon would be modified by 
the tensor interaction due to one pion exchange and more 
recent discussions may be found in e.g. the review pH 
and Ref. j^. A second effect is that the quasiparticle 
interaction contains explicit tensor contributions. Fol- 
lowing Refs. y. L|. it has generally been assumed that 
these have an exchange-tensor structure similar to that 
of the one-pion exchange interaction. However, Schwenk 
and Friman 7] have pointed out recently that the one- 
pion exchange interaction, when acting in second order, 
can give rise to contributions to the effective interaction 
which have a different structure. In their paper they 
evaluated these induced interaction contributions to the 
quasiparticle interaction and found that the exchange- 
tensor term in the quasiparticle interaction is much re- 
duced, and that terms of a different structure can be 
of a comparable magnitude to the exchange-tensor ones. 
A third effect of the tensor interaction is that there are 
multipair contributions to the magnetic susceptibility. 



2 



In this paper we begin by deriving a general expression 
for the magnetic susceptibility in Sec. 11111 taking into ac- 
count tensor contributions to the magnetic moment, and 
tensor contributions to the effective interaction which are 
completely general for an interaction which conserves 
the total angular momentum and the total spin of the 
quasiparticle-quasihole pair. This represents a general- 
ization of the earlier calculation of the magnetic suscep- 
tibility by Haensel and Jerzak |23j, who took into ac- 
count the effect of the exchange-tensor contribution to 
the quasiparticle interaction. Another issue that we ad- 
dress is how to parametrize tensor contributions to the 
quasiparticle interaction (Sec. IIV() . The scheme usually 
employed in the past suffers from the disadvantage that 
it is generally necessary to take into account high angu- 
lar harmonic contributions to the interaction. We present 
an alternative parametrization for which higher harmonic 
terms play little role. In Sec. we evaluate the matrix 
elements needed for the calculation of the susceptibility 
and estimate the magnitudes of the different contribu- 
tions. In Sec. I VII we extend the result to symmetric nu- 
clear matter and to responses with other spin and isospin 
properties, for example that to weak interactions. Sec. 
IVIII contains concluding remarks. 

II. BASIC FORMALISM 

The reason for Landau's theory of normal Fermi liq- 
uids being particularly simple for systems with inter- 
particle interactions that are central is that the ener- 
gies of long- wavelength, low-lying states differing little 
from the ground state may be described solely in terms 
of the quasiparticle distribution. The importance of con- 
servation laws in determining the low-frequency, long- 
wavelength behaviour of Fermi systems has been stressed 
by Leggett [24[ and a discussion in terms of Landau 
theory is given in Ref. [25|. Expressed in other terms, 
the only excitations of importance are ones with a sin- 
gle added quasiparticle and a single quasihole. When 
non-central interactions are present, the energy must 
include contributions from states with more than one 
quasiparticle-quasihole pair. To facilitate satisfying the 
conservation law for particle number, it is convenient to 
work with the thermodynamic potential (H— p,N) which, 
for brevity, we shall refer to simply as the energy. Here 
H is the Hamiltonian operator, N the particle number 
operator and fl the chemical potential. We write the 
change 8E in the energy when the system is excited from 
its ground state as the sum of a quasiparticle (Landau) 
contribution, 5El, and a multipair contribution, SEm'- 

SE = SE L + SEm- (2) 

The quasiparticle contribution to the energy may be 
expressed as the sum of two terms, the intrinsic kinetic 
and mutual interaction energies of the quasiparticles, and 
the energy of interaction with the magnetic field: 

5E L = 6E£ t + 6E?. (3) 



For simplicity, let us begin by considering a system of 
one species of fermions with spin 1/2. The change in 
the kinetic and interaction energies of the system when 
the quasiparticle distribution n pCT changes by an amount 
Sn p(T is give n to second order by the standard expres- 
sion [2(| |23 , which amounts to the statement that the 
quasiparticle contribution to the energy is given by 

(4) 

where e p is the quasiparticle energy. Note that Sn p (cr) 
is a matrix in spin space. The quantity / pp ' CTtr ' is the 
Landau quasiparticle interaction, which we write in the 
form 

/p P w = / P p' + .gpp'C • cr' + /p pVfflj (5) 

where / pp < is the spin-averaged quasiparticle interaction, 
3pp' is the spin-exchange contribution and / pp / CT(T / is the 
tensor or, more generally, the non-central contribution to 
the interaction. We shall discuss the form of the tensor 
interaction, / pp / CTCr /, in Sec. in the light of the new 
results of Schwenk and Friman 7] . To begin, we derive an 
expression for the magnetic susceptibility, for a comletely 
general interaction that conserves the total spin and the 
total angular momentum of the quasiparticle-quasihole 
pair, and subsequently we shall discuss its specific form 
in detail. 

Non-central forces can alter the effective charges of the 
quasiparticles. When calculating the energy of interac- 
tion of the quasiparticles with an external field it is there- 
fore important to allow for the fact that the effective cou- 
plings (charges) of quasiparticles may have components 
which are not scalars under rotations of the momentum of 
the quasiparticle. For definiteness, we consider the case 
of response to a magnetic field. The change in quasipar- 
ticle energy due to the application of a magnetic field, H, 
is given by 

Hij<TjHi6n p {<T) ^^ 3 , (6) 
where /xy, the magnetic moment matrix, is 

IMd = M% + i^— 2 g-J (7) 

where /1 and /it are parameters. Let us now consider the 
change in the energy when the Fermi surface is distorted. 
The distortion is specified by the function 

Kp)<*/3 = u (p) ' a aS3 (8) 

which corresponds to the change in the Fermi momentum 
as a function of direction. At zero temperature, the cor- 
responding change in the distribution function 5n p may 
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be expanded in powers of v and to second order one has 
(Sn p ) a/ 3 = S(p - p F )v(p) al 3 

- ls'(p-p F )^2^(p) a ^(p) ll3 . (9) 



It is advantageous to use spherical tensors so, following 
Refs. HE!, we write 



Kp) = E (-)"«"($)*" 

and expand u in spherical harmonics 



(10) 



(11) 



We shall work with eigenstates of the total angular 
momentum of the quasiparticle-quasihole pair. The cor- 
responding amplitudes are constructed by the transfor- 
mation 



Cfj = 52(-niml-n\JM)ujt m , 



(12) 



where (Zml — (i\JM) are Clebsch-Gordan coefficients. 

The part of the quasiparticle energy not involving the 
external field, Eq. (0J, has been calculated in Refs. [5ll23j 
to be 



5E\ nt 



Pf 



{2nh) 3 m 



7 E C%Cft{l'J\A\lJ), (13) 



Rewriting this in terms of the coefficients Cfj and adding 
the intrinsic contribution to the energy leads to the fol- 
lowing expression for the quasiparticle energy: 



SE L = 



Pf 



(2Trh) 3 m 
2m 



E C%CM*{l'J\A\U) 



Pf 



ll'JM 

AttH 



1 y/2 1 



At first sight, one may be surprised that there are no 
terms with J = 3. However, this follows immediately 
from the observation that the magnetic moment operator 
[lijOj transforms as an axial vector under simultaneous 
rotations in momentum space and spin space, and con- 
sequently it contains only components corresponding to 
total angular momentum J = 1. 

By minimizing the quasiparticle contribution to the 
energy with respect to the coefficients Cfj we find that 
all Cfj are zero except the following two, 

Cq ! = V^H Z — (n(21\A\21) + ¥±{01\A\21) 
PF V V2 

x , (17) 



(21|A|21)(01|A|01) - |(01|A|21)f 



and 



C° 21 = V^H Z — (-M<01|il|21)-^(01|A|01) 
Pf V V2 

(21| J 4|21)(01| J 4|01)- |(01|A|21)| 2 ' {lS) 



ll'JM 



where the matrix elements in the new basis are given by 
(l'J'\A\U) = E (J'M'\l'm'l-n') 

(l'm'fi'\A\lmfj)(lml - n\JM), (14) 

{V J'\A\IJ) are the normalized matrix elements, and 

A=l + T, (15) 

the unit coming from the first term in Eq. Q, and 
J- = ^(0)/ P p' CTCT ' comes from the interaction term. For a 
quasiparticle interaction containing only exchange-tensor 
terms, the matrix elements (U\A\l'J) have been calcu- 
lated for the conventional parametrization for the ex- 
change tensor interaction in Ref . [j| . 

If we take the magnetic field to be in the ^-direction, 
TL = Ti Zl we find that the energy due to the application 
of this external field is 

8E? = 2 4^n z x 



{2iThf 



M^oo + "^"20 - g Ki + U 2-l) 



.(16) 



III. STATIC RESPONSE FUNCTIONS 

The magnetic susceptibility is defined as the derivative 
of the magnetization with respect to the magnetic field 
at zero field, 



dM,, 



X 



(19) 



We write the magnetization, M, as the sum of 
two terms, one contribution from single quasiparticle- 
quasihole pairs, Ml, and another contribution from the 
excitations of multipair states, which we denote by Mm- 



M = M L + Mj 



M ■ 



(20) 



The first term, the Landau term, may be calculated by 
the same methods as we used in the previous section, and 
in terms of the CM it is given by: 



Mr. 



2p 



{2nh) 3 



\71 



(21) 
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By using Eqs. JTJJ and JTHJ for the Cfj, and the rela- 



tion between the susceptibility and the magnetization, 
we obtain the following expression for the susceptibility: 



X 



N(0) 



/4 



(21|A|21) 



(21|A|21)(01|A|01) 
(01|A|21) 



(01|A|21) 



(21|A|21)(01|A|01) 
(01L4|01> 



(01|A|21> 2 



2 (21|^|21)(01|^|01) 

XM, 



<01|A|21> 



(22) 



where xm is the contribution from transitions to mul- 
tipair states. Here we have implicitly made the usual 
choice of phase for the states, in which case the matrix 
elements of A are real. If we take only the first term and 
set /x = /x in Eq. 121'P , this result reduces to that of Ref. 
pflj . Quantitative estimates of the different effects will 
be described in Sec. E] after calculation of the relevant 
matrix elements, but first we shall discuss the form of 
the tensor interaction. 



IV. THE TENSOR INTERACTION 

The explicit form of the tensor interaction including all 
terms allowed by invariance under simultaneous rotations 
in coordinate space and spin space is 

= Vp'S(q) + k pp ,s(P) + v^(q, p)- (23) 

Here 



Note that the choice of parametrization of the exchange- 
tensor term in Eq. (|23|) is different than the one tradi- 
tionally used. The conventional way to parametrize this 
term is to write \A 



,S(q) 



<r_i 

2 "PP' 



5(q). 



(28) 



The functions h and h contain the same physical infor- 
mation, but have different properties as far as their ex- 
pansions in terms of Legendre polynomials are concerned. 
We will now compare some aspects of the different choices 
of parametrization. 

With the conventional parametrization, h pp i, the ex- 
plicit q 2 factor ensures that the tensor interaction van- 
ishes identically for q — > 0, provided that h pp i is finite 
for p — > p'. This is in agreement with the form of the 
one-pion-exchange contribution which is given by 



T 

pp'ercr' 



3m 2 q 2 



S(q), 



(29) 



where / is the pion-nucleon coupling constant, and 
f 2 / (Antic) ~ 0.08. More generally, if the tensor inter- 
action is analytic for small q, one would expect h pp > to 
tend to zero for p —> p', since otherwise the tensor inter- 
action would depend on the way in which p approaches 
p'. Such nonanalytic contributions to the quasiparti- 
cle interaction do occur for systems with central inter- 
actions [28L l29l| , and there they are due to exchange of 
quasiparticle-quasihole pairs and collective modes with 
long wavelengths. The work of Schwenk and Friman 
indicates that such contributions do not occur in h pp i, 
although they do in k pp i (see Sec. IVBfl . To ensure that 
hpp' ~~ * for p — » p', the hi must satisfy the condition 



S(q) — 3<r • q er' • q <x • <x' 



(24) 



is referred to as the exchange-tensor operator, q = p 
p' being the momentum transfer. The operator in the 
second term, which is referred to as the center of mass 
(cm) tensor, is given by 



S(P) = 3<r • P a-' ■ P - tr • a'. 



(25) 



It has the same form as the exchange-tensor, but with q 
replaced by the total momentum P = p + p'. The final 
term, referred to as the cross-vector term, is given by 

A(q,P) = (<r x tr') ■ (q x P) 

= a ■ P a' ■ q- a ■ q a 1 P. (26) 

The functions h pp i, k pp > and / pp ' may be expanded in 
Legendre polynomials of cos#, where 9 is the angle be- 
tween p and p', e.g. for h pp >, 



h pp > = ^ hi Pi (cos ( 



(27) 



(30) 



1=0 



If the series is truncated after a finite number of terms, 
this will inevitably lead to the sum rule being violated. 
However, this should not be very important as long as 
long- wavelength processes do not play a dominant role. 

The conventional parametrization suffers from the dis- 
advantage that, when expanding h in terms of Legen- 
dre polynomials l|27() , it is generally necessary to include 
many terms in the sum in order to obtain a good rep- 
resentation of the interaction. Consider, for example, 
the one pion exchange contribution. Using the fact that 
COS0 = 1- q 2 /(2p 2 



Op) one finds 



h 



OPE 



f 2 r 2pF da 2 1 



(31) 



i^ m l) Jo Pf I 2 + m l 



1=0 



where Qi is a Legendre function of the second kind. In 
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TABLE I: Tensor parameters for neutron matter with Fermi 
wave vector A?f = 1.7 F . Values for the tensor parameters 
by Schwenk pjoll . obtained taking into account the tensor in- 
teraction to first order by the renormalization group method 
described in Refs. 0,1^1 are presented in the first and third 
columns, for the two choices of parametrization. For compar- 
ison we present the result for the one-pion exchange potential 
(Eq. 13210 in column two and four. All tensor parameters are 
calculated for m*/m = 0.8345. 



1 


H, 


ff OPE 


Hi 


£OPE 





0.529 


0.403 


0.665 


0.751 


1 


0.150 


-0.088 


1.200 


1.037 


2 


-0.0959 


-0.079 


1.150 


0.935 


3 


-0.141 


-0.064 


0.933 


0.753 


4 


-0.124 


-0.049 


0.703 


0.575 


5 


-0.0944 


-0.036 


0.509 


0.427 


6 


-0.0697 


-0.026 


0.359 


0.310 


7 


-0.0494 


-0.019 


0.249 


0.223 


8 


-0.0335 


-0.013 


0.171 


0.158 


9 


-0.0225 


-0.009 


0.116 


0.111 


10 


-0.0151 


-0.007 


0.0774 


0.078 



Table 1 we show results of calculations of the exchange 
tensor interaction, Hi = N(Q)hi, by Schwenk [3(J, for 
fc F = 1.7 F _1 . Note that these were calculated taking 
the tensor interaction into account only in first order. 
One sees that as I increases from zero, Hi first increases 
and then decreases. It is therefore a poor approxima- 
tion to truncate the expansion after the first few terms. 
The reason for this is that typical Fermi wave numbers 
at nuclear densities are of order 1.5 F _1 , which is large 
compared with the pion mass, m^/hc ~ 0.7 F _1 . In the 
limit of the zero pion mass, all the hi diverge, the leading 
contribution being (21 + l)f 2 /(3m 2 ) \n(2pp / m^) . This is 
due to the fact that h ~ q~ 2 , which leads to a logarithmic 
integral for hi. 

The situation is quite different for the parametrization 
(|27p . For the one pion exchange potential (|29[1 one finds 



f 2 



(3m*) 



2p| 



l(2l + l)Q l [l + m 2 j2p F ] .(32) 



This shows that for small m^/pF, the leading term is 

JL 2 A± X ( m * 

1 ~ 3ml 2 



— hy(2p F /m w ). (33) 

PF 



In the limit of zero pion mass, h pp i becomes a constant, 
and therefore the only nonvanishing coefficient in the ex- 
pansion in Legendre polynomial is ho . Note that for both 
hi and hi, the limits m^/p-p — > and I — > oo do not com- 
mute. 

Let us now compare the expansions of hi and hi . Using 
the relation between h pp i and h pp i, Eq. I|28|) and that 
1 — cos# = q 2 /(2p F ), we find the following expression for 



hi = 2(21 + 1) f + ^ C ° S ^ (1 - cos60Mcos6>)P;(cos 6>){34) 



. I + 1 - I - 
hi = 2 [ hi - — i n hi+i - — -hi-i 



21 + 3 



21 - 1 



(35) 



Thus, given the hi's, we can obtain the /i;'s. Next we 
show how to write the hi's, given by 

- (21 + 1) f 1 d(cos0) h(cos0) 
h = — 7. — / — 7. — ; — - p i( cosd ) (36) 



2 y_! 2 l-cos( 



in terms of the hi's. Expanding the function h(cos8) in 
Legendre polynomials, we get 



hi = 



(21 + 1) / ,1 d(cos 



• 5^/1!/ fl/ (COS fl)- 



Pi(cos6) 



cos 9 



(37) 



Since the series on the right hand side is not uniformly 
convergent for cos 8 in the closed interval [—1, 1], we can- 
not invert the order of the sum and the integral. However, 
because of the sum rule (|30|l . h(6 = Q) vanishes, there is 
no divergence, and we can rewrite Eq. (|36|) as 



hi _ (21 + 1) ^ rfCcosfl) (fo(cosfl) - h (0 = 0)) Pi{coze y 



1 — cos ( 



Making use of the sum rule, we find that 
~ (21 + 1) f 1 d(cos( 



(38) 



■53/»,,[^/(oose)-l]- 



Pi(cos8) 



— cos ( 



(39) 



This expression may be evaluated using a number of stan- 
dard results. First 



1 d(cos9) [P v (cos 8) - P (z)} 



Pi (cos ( 



P v (z)Qi(z)-P (z)Qi(z), l'<l. (40) 
The Qi(z) may be expressed as 

Q l ( z ) = l Pl ( z )ln\±^-Wi- 1 , (41) 
2 1 — z 

where 



' 1 

Wi-i = ^-P fc -i(z)P ; _ fe (z) and (42) 
fe=l 

W-i = 0. (43) 
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The Legendre polynomial Pi(z) may be written in the 
following form: 

flW = 2-'g(-l)'f ^ ( 2/ ; 2fc )^ 2fc - (44) 

Letting z — > 1 in the integral, we find that the logarithmic 
divergence vanishes and the final result is 

r (21 + 1) ^ i 

l' =1+1 k=l + l 

V. EVALUATION OF THE MATRIX 
ELEMENTS 

We return to the calculation of the matrix elements of 
A = 1 + J 7 , Eq. (| 1 5p . which we write as 

(ZJ|A|Z'J) = + (U\F\l'J). (46) 



We further write (iJ|J r |i'J) as the sum of a number of 
contributions: 

(U\F\l'J) = (U\F c \l' J) + (U\F c *\l' J) 

+ (IJ\F CVCI \1'J) + (IJ\F CV \1'J). (47) 

Here T c is the central part, JF CX is the exchange tensor 
part, J- cm is the center of mass tensor part and J- cv is 
the cross vector part. 



A. Central and exchange tensor contributions 

In Ref. the central and exchange tensor part in the 
conventional parametrization were calculated to be 

(U\F c \l'J) = <^'2lTT' (48) 

and 



(U\F*\l' J) = 5 U 



21 



(-1) J 15 



1 1 2 



v 

3[(2/ + l)(2/' + l)] 1 /2 ( i i j 



2J+1 



(21 - 1)(2Z + 1) 

I V 2 \ ( 1 V J 
/ 1 I 1 2 



Hi 







i+1 



Hi 



2(1 + 1) 



(2Z + 3)(2/ + 1; 




21 + 1 



2V 



1 



1/2 



Hv 



2V + 1 
21 + 1 



1/2 



iJ J -3[(2Z + l)(2Z' + l)] 1 / 2 ^^ 




1 V J 
1 I I" 



(49) 



Here 



h h I 
m\ m 2 m 



are Wigner 3j-symbols and 



are Wigner 6j-symbols. 



I h h h 

{ h ^2 h 

We can transform the matrix elements above to the 
new parametrization. In the expression 1221) we can see 
that the only matrix elements that are necessary for the 
calculation of the susceptibility are (01|A|01), (01|A|21) 
and (21|A|21). The contribution to these matrix elements 
from non-interacting quasiparticles and from the central 
part of the quasiparticle interaction is given by 

(l'j\l+r\U)=5 w (l+2^l) (50) 

and the exchange tensor contributions by [j| 
(Ol^Ol) = 0, 

(01|^ ox |21) = -V2 (h - ^Hi + ^H 2 ^j and 



(21|^ cx |21) = -Lh 1 + Ih 2 -^-H 3 . (51) 
15 5 35 



Let us now express this result in terms of Hi. As an 
example we consider the combination Hq — 2Hi/3 + H 2 /5. 
This may be written in the form 



Ho 



2 -Hi 
3 



l -H. 

5 



+ld(COS0) ^(cos^)[P o (cos ( 
+1, 



2 

l d{cm 9) 



H(cost 



[Po(cos( 



2Pi(cos6 
2Pi(cos( 



P 2 (cos t 
-P 2 (cost 



2(1 



(52) 



An interesting property of the fraction in Eq. (|52|) is that 
it is finite for 6^0, and therefore the integral converges 
if H(cos9) is finite in this limit, even though individual 
terms diverge. One finds similar results for the other 
combinations of the Hi that occur in Eqs. J5TJ. The 
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corresponding result written in terms of the Hi are 
(01|JF ex |01) = 0, 

(01\F>*\21) = -^=(H -H 1 ) and 

(2i|^»|21) = ^-i(^ + ^i). (53) 

It remains to calculate the center of mass tensor and the 
cross vector contributions. 



B. Center of mass tensor 

The calculation of the center of mass tensor terms fol- 
lows closely that of the exchange tensor. Let us study 
the operator JF cm given by 

= N(0)k pp ,S(P) = KiPi(cos9)S(P). (54) 



will be finite for 8 — > ir, provided that fc pp ' is finite in 
this limit. Using the result of Schwenk and Friman 
that fcpp/ goes to a constant when 8 — > ir, we can see that 
hi, which is given by 



h = 



21 + 1 



1 + cos 



-Pi(cos8)d(cos8), (58) 



will, in fact, have a logarithmic divergence. To deal with 
this divergence, we replace the limit —1 by A, and write 
ki as the sum of a regular part and a part that contains 
the logarithmic divergence, 

fei(A) = fc[ cg + (-l)'(2i + 1) ln(l + A)C(A) (59) 

where C(A) = fc pp /(A)/(2p|) tends to a constant when 
A — » —1 and the factor (—1)' comes from P;(— 1) = 
(—1)'. By inserting the expression above in the matrix 
elements of JF cm , Eq. ljBT)l. we obtain 



If we reverse one of the momenta, we find that 



-77cm 
p-p' 



N(0)k p ^ p ,S(q) 
J2(-^ l KiPi(cos6)S(q). 



(55) 



Thus we have a similar expression to the exchange tensor 
operator, with a difference of a factor of (—1)'. For the 
state vectors \l'J'), reversing the momenta, gives 



-p'a'\l'J') = (-l) l '{p'a'\l'j'). 



(56) 



The matrix elements of J- cm is thus the same as those for 
J-" ox , apart from an overall factor of (—1)' and that Hi 
should be replaced by (— 

However, when calculating the matrix elements of the 
center of mass tensor, it is necessary to be careful for 
8 — ► 7r, since the calculation of fc pp ' by Schwenk and 
Friman Q showed that k pp i does not vanish in this limit, 
but rather tends to a constant when second order tensor 
contributions are included. However, we shall show, that 
for the matrix elements we consider, this will not pose a 
problem. 

By writing fc pp / as {P 2 /pp)k PP < we can use the cal- 
culation for the exchange tensor given in Ref. [Bj], ex- 
cept that for the odd I' states we have to include a fac- 
tor (—1)' , and Hi should be replaced by (—l) l Ki, with 
Ki = N(Q)ki. For the particular matrix elements that 
we are interested in we find that 

(01|^ cm |21) =-V2^k + ^K 1 + ^K 2 ^j and 
(21|;r<-|21) = ^K, + \k 2 + ^K 3 . (57) 

Just as in the example (|52|l . we find that the matrix ele- 
ments above give rise to combinations of Legendre poly- 
nomials which contain a factor 1 + cos 8, so the integral 



(01|.f cm |21) 



V2 (i^ og + C(A) ln(l + A) 
2C(A) ln(l + A)) 



3 1 



+ -K r 2 cs + C(A) ln(l + A) 
o 



-V2^ + ^re g+ l^re g | ((i()) 



and 



(21|jr cm |21) 



^[i^ cg -3C(A)ln(l + A)] 

-K r 2 e£ + 2C(A) ln(l + A)l 
5 

l[^ eS -7C(A)ln(l + A)] 
7 ~ 9 ~ 3 - 

15 Kl + 5^ 2 + 35^ 3 • 



(61) 



This shows that in the combinations of parameters we 
need, the divergent contribution vanishes and therefore 
allows us to simply replace Hi by (— 1) Ki in the expres- 
sion £3) for the (U\T cyL \V J) to obtain \lJ\T cva \l' J) in 
terms of the Ki's: 

(OlI^lOl) =0, 

<01|^ cra |21) = -^(^ + ^ 1 ) and 
(21\^\21) = ^-l(K -K 1 ). (62) 



C. The cross vector term 

Finally, we turn to the cross vector term, which is pro- 
portional to 



A(q,P) = (a x a') ■ (q x P) 
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= O x a') ■ 



P x P 
|p x P' 



er • per' • p — <r ■ p'er' ■ p 

|pxp'| 



(63) 



The calculations of Ref . Q indicate that this term will be 
well behaved in the limits where sin# — > 0, since L pp i — > 
for both 9 = and 9 = it. This means that there will 
be sum rules like the one (|30|) for both 9 = and 9 = it: 

^ZjP;(cos0 = 1) =0, and 

j2 L i p i( cos ° = - 1 ) = X^ -1 )'^ 00861 = ^ = 

which gives two sum rules, one for the odd coefficients 
and one for even ones: 



i=a 



0, and 



E L21+1 
1=0 



= 0. 



(64) 



To simplify the calculation of the matrix elements, we 
define L pp > = L pp i / sm9, and calculate the matrix ele- 
ment (I J I T CV \V J), which later will be transformed back 
to the original notation. The cross vector contribution 
to the quasiparticle interaction is 



AF CV = N(0)~Tr, 

47T 4 L 



dndfL'u ■ ctT^m 1 ■ er'.(65) 



We use the notation and method described in Sec. ITT1 
which gives for the matrix elements in the \lml — fj) basis 



(Iml-nlJ^ll'm'l-n') = 



I 1 I" \ V 1 I" 

A 



/ 

m 



-fj, m 



I" 



E 

' V 



(-) r 



"l. 



1 1" 

-H' m" 



I 1 I" 




r 



v 1 l" 
000 



1 1 i" 

m —fx' m" 



(66) 



V 1 I" 
m' —fx m" 



Changing to the \IJ) basis, we find the matrix element 
to be 



(U\F cv \l'J) = 

[(2/ + i)(2r + i)] 1 /2 / 1 1 jVi' 1 j 
I A 



2J+ 1 
[(2I + 1)(2I'-F 

' I 11" 




L., 



El 



i/r 1 r 
)[ 



1 /' j 
1 i i" 



(67) 



For I' = I + 2, the triangle conditions imply that the 3-j 
symbols are non-zero only if J = I + 1 and Z" = I + 1. 
When this happens the two terms are equal in magnutidc 
but of opposite sign. Therefore all (U\!F cv \(l + 2) J) are 
zero. The only non-zero matrix elements are when I = V . 
Then J can be either I + 1 or Z — 1: 



(ZJI^IZJ) = 
2(Z + 1) 2 
(2Z + 3) 2 (2M 
Z 

(2Z + 1)(2Z- 
(ZJ|^ V |ZJ) = 



2) 



-L 



1 



H)(2J + 1)-1_ 
for J = Z + 1, and (68) 



I 



(2Z-1) 2 



1 



Z(2Z + 1) 



2(Z + 1) 2 



(2Z + l)(2Z + 2)(2Z + 3) 



L1-1 

Li+i, ior J = 1-1 (69) 



The matrix elements we need thus have the following 
contributions from the cross vector: 

(01|^ cv |01) = 0, 
(21|^ cv |01) = and 

(21|^ cv |21) = iLx-lia- (70) 

The only matrix element we need for the cross vector 
term is thus when I = V = 2 and J = 1, and we transform 
this to our original parametrization by writing: 



£d(oos0)|x; 



which is the integral 



21|JT CV |21) 
L k P k (cos( 



sm( 



(71) 

[Pi (cos 6) -P 3 (cos 6)] 



f dxJ2 L nPn(x)Vl -X 2 X = ~ ^i„4- (72) 
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TABLE II: Landau parameters for neutron matter calcu- 
lated by Schwenk (Tf.The Landau parameters are calculated 
for Fermi wave vector k& = 1.7 F _1 and the effective mass is 
m*/m = 0.8345. 



1 


G, 


Ht 


Ki 


Li 





0.842 


0.070 


-0.258 


-0.060 


1 


0.412 


0.163 


0.146 


-0.089 


2 


0.219 


-0.301 


0.124 


-0.034 


3 


0.109 


-0.002 


0.048 


0.025 


4 


0.051 


-0.150 


0.015 


0.043 



pure neutrons, we first use the Landau parameters ob- 
tained by Schwenk |3(j listed Table [I] These take into 
account the tensor interaction only to first order, and 
consequently the only tensor term is the exchange ten- 
sor. We find 



We calculate the integral /„ to be 

( [n/2] 



In = < 



E 

k=Q 




7r(-l) fc (2n - 2Jfe)!(n - 2k + 1) 
k\( n - jfe)!(s=»±i )i(»=»±3jr 



, n = 2l + l 



n = 2l 



D. Matrix elements and the susceptibility 



Combining Eqs. (JSDJl , (£23) , © , G2J> and G2J) above, we 
find 



(01|A|01) = 1 + G , 

(01|A|21) - [£T - If! + K + Ki] 

<21|A|21) = l + i( G2 + ^±^ 

(ff + ^o + Hx-Kx)-\^L n I n , 



and 



(74) 

with I n given by Eq. I|73|) . 

From these results one may calculate the magnetic sus- 
ceptibility from Eq. I|22|). The effects of the tensor in- 
teraction are generally small, so to estimate the various 
contributions to the susceptibility we expand the suscep- 
tibility to second order in the tensor interaction and (it 
and find 



N(0) 



1 



1 



<01|A|21) 



a4 



1 + G ' " (1 + G ) 2 1 + G 2 /5 

(l + G )(l + G 2 /5) 

XM- (75) 



2 l + G 2 /5 



It is not easy to obtain a reliable estimate of the various 
terms, since different calculations of the tensor Landau 
parameters give different result. To exemplify this for 



(01|A|21) ps -0.13 
(21|A|21) ps 0.86. 



and 



(76) 



We can now compare this with what we get if we use 
the Landau parameters of Ref . , listed in Table [H] 
These were calculated including the tensor interaction to 
second order. We find that the exchange tensor, center 
of mass tensor and the cross vector terms are all quite 
small, so the matrix elements are only slightly modified 
from the value one obtains using only the central part 
of the interaction. In addition to being small, the tensor 
contributions have a different sign from the first order 
result: 



(73) 



(01|A|21) re 0.07 and 
(21|A|21) ps 1.06 with 
(21|1 + .F C |21) PS 1.04. 



(77) 



Since the calculation of tensor parameters in Ref. Q is 
only to second order, and the second order effects almost 
completely cancel the first order ones, higher order effects 
could be significant and should be calculated. 

From the two sets of results above we can draw some 
conclusions regarding the effect of the tensor interaction 
on the susceptibility. The second term in Eq. I|75l) is 
clearly small compared to the first term, of the order of 
1%, and may be neglected. The third and fourth terms, 
could be important, depending on the value of the tensor 
magnetic moment. The corrections to the bare nucleon 
moment arise from two sources. The first is configuration 
mixing, which is due to the fact that a quasiparticle in 
the medium consists of a bare nucleon and a superposi- 
tion of more complicated states involving a nucleon plus 
a number of particle-hole pairs. The second source is 
exchange currents, which are due to the fact that the ex- 
ternal field can interact with other degrees of freedom 
than nucleons, for example, intermediate mesons in a 
nucleon-nucleon interaction. For definiteness, we shall 
use the results given by Arima et al. in Ref. [2l| for 
the Fermi gas model. For neutrons in symmetric nuclear 
matter, they find fi ~ —1.62 /in and fj,T ~ 0.08 /xn, 
where /^n = eh/2M is the nuclear magneton. Calcula- 
tions of the configuration mixing contribution to weak 
interaction matrix elements made by Cowell and Pand- 
haripande |32| for asymmetric nuclear matter indicate 
that the renormalization of the magnetic moment does 
not depend strongly on proton fraction, and amounts to 
a suppression by about 10%. We would thus suspect that 
renormalization of the magnetic moment would be rela- 
tively insensitive to the proton fraction, and would be of 
a similar order of magnitude. However, we stress that it 
is important to make more detailed estimates of magnetic 
moments in the nuclear medium. 
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The small magnitude of the tensor magnetic moment 
means that we may neglect also the third and fourth 
term in Eq. (|75|l , both being less than one per cent of 
the first term. Therefore, we conclude that, the most 
important effects of the tensor force on the susceptibility 
is most likely the renormalization of the isotropic part of 
the magnetic moment and to the presence of transitions 
to multipair states. 

No reliable calculation of xm exists, but by using the 
sum rule argument presented in Ref. |19| we can put a 
lower bound on the multipair contribution to the suscep- 
tibility: 

XmM)>^?M, ( 9 -0). (78) 

Here S(q) is the static structure function, n is the den- 
sity and Q is the mean excitation energy. In Ref. [l9l | 
we took the values of S(q) = 0.19 and u> = 63 MeV from 
the calculations of Ref. [3i| and estimated the multi- 
pair contribution to the susceptibility to be more than 
60% of the susceptibility calculated by Fantoni et al 
(x = 0.38xf)- Here xf is the susceptibility of a free 
Fermi gas of neutrons, xf = 3n/2ep- New calculations 
of the static structure function and the mean excitation 
energy by Cowell and Pandharipande |34j give a similar 
result, S(q = 0) = 0.187 and Q = 64 MeV. In order 
to compare with symmetric nuclear matter in next sec- 
tion we calculate xm/xf, and we find it to be ~ 0.23 for 
neutron matter, if we use the calculations of Ref. [34|. 

Finally, we note that the calculation for the total sus- 
ceptibility of Fantoni et al. ^3 gives a result very close to 
the value one get by calculating G from Brueckner the- 
ory and using Eq. JJ| for the susceptibility, which does 
not take into account the renormalization of the mag- 
netic moment. This would appear to indicate that the 
tensor correlations redistribute spectral weight between 
single pair and multipair excitations, leaving the total 
susceptibility unchanged. A similar effect is familiar for 
the compressibility of electrons in metals when one takes 
into account the electron-phonon interaction [24[ . 



VI. SYMMETRIC NUCLEAR MATTER 

Symmetric nuclear matter can be treated by a straight- 
forward generalization of the discussion above by intro- 
ducing the isospin degree of freedom in addition to the 
spin one. Thus the generalization of Eq. JSJ becomes 

/ PP w = /p P ' + /pp/T • t' + g pp >cr ■ ct' + g pp ,er ■ a'r ■ t' 

+ /ppW+/ppWT-T', (79) 

where r is the isospin operator. Since the magnetic 
moments of neutrons and protons are different, a mag- 
netic field couples to both isoscalar and isovector com- 
ponents. Because the difference between the magnetic 
moments of a neutron and a proton is much larger than 



TABLE III: Landau parameters for symmetric nuclear matter 
calculated by Q. The Landau parameters are calculated for 
Fermi wave vector hp = 1.35 F _1 and the effective mass is 
taken to be m* /m = 1. 



1 


Gi 


G\ 


H t 


Hi 





0.447 


1.291 


0.65 


-0.255 


1 


0.760 


0.070 


0.975 


-0.359 


2 


0.276 


0.090 


0.829 


-0.315 


3 


0.078 


0.069 


0.617 


-0.233 



the sum, the magnetic response will be dominated by 
the isovector term. The isoscalar spin response is given 
by Eq. J22J), but the density of states to be used is 
iV(0) = 2m* py /Tr 2 h 3 ', the factor of two reflecting the two 
nucleon species. The isovector spin response is given by 
a similar expression, but with G, H, K and L replaced by 
G',H',K' and V . Since there is no calculation of the 
new tensor terms for symmetric nuclear matter, we esti- 
mate the matrix elements by using the exchange tensor 
Landau parameters for symmetric nuclear matter calcu- 
lated in Ref. p|, which are listed in Table ITTT1 This gives 
the contribution to the matrix elements to be 

(01|A'|21) w 0.11 and 

(21|A'|21) sa 1.08. (80) 

The calculations of Ref. [2lJ give for the anisotropic 
contribution to the isovector magnetic moment the value 
/it ~ 0.08 fi n and therefore the third and fourth terms 
in the analogue of Eq. (|75|l for the isovector magnetic re- 
sponse may be neglected. The second term is also unim- 
portant, since it is on the order of 1% of the first term. 
We now use sum-rule arguments |l9j to put a bound 
on the fifth term, which comes from multipair excita- 
tions. Equation l|78[) applies to a Fermi system with an 
arbitary number of components, and recent calculations 
by Cowell and Pandharipande [34| of the static struc- 
ture function and mean excitation energy for the isovec- 
tor spin response. Their calculations give S(q) = 0.155 
and u) « 253 MeV in the limit q — 0. Comparing the 
contributions from multipair excitations with the suscep- 
tibility of a free Fermi gas we find 

Xm(9,0) > 4e F 5'(q) ^ 
Xf ~ 3d) 

where xf = 3n/2ep is the susceptibility of a free Fermi 
gas consisting of 2 species of spin- 1/2 particles and e F 
is the Fermi energy. At nuclear matter density n = 
0.16 fm~ 3 , cf ~ 37 MeV, gives that xm is at least ~ 3 % 
of xf, which is much smaller than the corresponding 
bound for pure neutrons. 

We are not aware of recent calculations of the total sus- 
ceptibility for symmetric nuclear matter similar to the 
calculations of Ref. [l^, but we note that for neutron 
matter those calculations gave a susceptibility very close 



11 



to the result one obtains by calculating Go from Brueck- 
ner theory and using Eq. (JIJ to calculate the suscepti- 
bility. Therefore we compare xm with the analogue for 
symmetric nuclear matter of Eq. l(T]l. which we call xo ; 
and find 

XM^4 £F %) (1 + GT ) m 
Xo 3w m*/m 

where we have taken m*/m = 0.8, and G' from Ref. 
0, as listed in Tabic ITTT1 This result indicates that the 
multipair contributions, being at least 9 % of the total 
susceptibility, can have some importance for symmetric 
nuclear matter, even though the bound is much smaller 
than the corresponding bound for neutron matter. It 
would be valuable to have explicit calculations of the 
multipair contributions to the susceptibility since it is 
unclear to what extent the bound provides a realistic es- 
timate of the value of the quantity. 

VII. CONCLUDING REMARKS 

In this paper we have derived an expression for the 
magnetic susceptibility of a Fermi liquid for a general in- 
teraction which conserves the total angular momentum 
and the total spin of the quasiparticle-quasihole pair. 
Apart from the contribution from single quasiparticle- 
quasihole pairs, which may be calculated using Landau 
theory, the susceptibility also contains a contributions 
from excitation of multipair states. In addition to the 
exchange tensor term usually included, we have taken 
into account the new tensor terms found in Ref. Q. We 
have also introduced an alternative paramctrization of 
the exchange tensor interaction which has the advantage 
of reducing the importance of high angular harmonic con- 
tributions. 



For neutron matter we find that the most important 
effects of the tensor force are to renormalize the magnetic 
moment and to give a contribution from multipair states. 
For symmetric nuclear matter one can draw the same 
conclusion, the bound on the multipair contribution is, 
however, much smaller, but it is not clear how close the 
actual contribution to the susceptibility, from multipairs, 
is to this lower bound. 

The formalism may easily be extended to calculate the 
response to weak probes: the only difference being that 
the magnetic moments must be replaced by the corre- 
sponding matrix elements of the weak charges. 
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